arXiv:1503.00838v2 [cond-mat.mes-hall] 1 Feb 2016 


Majorana fermions in a superconducting Mobius strip 

Yuan Pang, 1 ’^ Jie Shen, 1 ’^] Fanming Qu, 1 ’^] Zhaozheng Lyu, 1 Junhua Wang, 1 
Junya Feng, 1 Jie Fan, 1 Guangtong Liu, 1 Zhongqing Ji, 1 Xiunian Jing, 1,2 
Changli Yang, 1,2 Qingfeng Sun, 2,3 X. C. Xie, 2,3 Liang Fu, 4 and Li Lu 1, 2, [£] 

1 Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, 

Beijing 100190, People’s Republic of China 
2 Collaborative Innovation Center of Quantum Matter, 

Beijing 100871, People’s Republic of China 
3 International Center for Quantum Materials, 

Peking University, Beijing 100871, People’s Republic of China 
^Department of Physics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 02139, USA 
(Dated: February 2, 2016) 



1 



Abstract 

Recently, much attention has been paid to search for Majorana fermions in solid-state systems. 
Among various proposals there is one based on radio-frequency superconducting quantum interfer¬ 
ence devices (rf-SQUIDs), in which the appearance of 4-7r-period energy-phase relations is regarded 
as smoking-gun evidence of Majorana fermion states. Here we report the observation of truncated 
47r-period (i.e., 27r-period but fully skewed) oscillatory patterns of contact resistance on rf-SQUIDs 
constructed on the surface of three-dimensional topological insulator B^Tea- The results reveal 
the existence of 1/2 fractional modes of Cooper pairs and the occurrence of parity switchings, both 
of which are necessary signatures accompanied with the formation of Majorana fermion states. 
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Majorana Fermions are a mysterious type of particles predicted eight decades ago 1 but 
are still nowhere to be found. Recently there is a hope to find Majorana quasiparticles 
in condensed matter system^ 2 " 4 . The pairing of helical electrons in topological insulators 
(TIs) and related materials are expected to form topological superconductivity resembling 
that of a spinless p -wave superconductor, and hosting Majorana fermion stated with which 
topological quantum computers could potentially be built. One of the signatures of Majorana 
fermion states, the appearance of a zero-bias conductance peak at the superconductor-normal 
metal (S-N) interface, has already been observed experimentally 6 . Also highly-sought is 
the appearance of 47r-period energy-phase relations (EPRs) in S-Tl-S type of Josephson 
junctions' 11 . To search for the 47r-period EPRs, a number of phase-sensitive experiments 
have been conducted on various Josephson devices constructed on 2D and 3D TIs 12 but 
clear evidence is still awaited. 

In 2008, Fu and Kane 5 proposed that Majorana fermion states should exist in Josephson 
junctions constructed on the surface of a 3D TI, as long as the phase difference across 
the junction is 7r. To maintain the i r phase difference, an ingenious way is to make a rf- 
SQUID by connecting the Josephson junction with a superconducting ring (Fig. 1) and 
threading half flux quantum into the ring. Such an experiment has been proposed in several 
different versions by the theorist^" 11 . For rf-SQUIDs constructed on the surface of 3D TIs, in 
particular, Wieder, Zhang and Kane predicted that the existence of Majorana fermion states 
in the junction would accompany with a 47r-period signature in the channel conductance of 
the junction, resulting in quantum phase transitions and conductance jumpings at every 
odd multiples of half flux quantum 11 . In such case, we anticipate that the local conductance 
spectrum should also oscillate, which should be detectable by using nano-probes placed on 
the surface of the 3D TI near the Josephson junction, as illustrated in Fig. 1. 

We have fabricated about twenty such devices by depositing Pb rings on exfoliated Bi 2 Te 3 
flakes of ~100 nm in thickness, and measured the contact resistance between normal-metal 
Pd nano-electrodes and Bi 2 Te 3 at positions A, B, C and D as marked in Fig. 1, by us¬ 
ing standard three-terminal measurement configuration and lock-in amplifier technique at 
cryogenic temperatures of dilution refrigerators. We found that the contact resistance at 
positions A and B oscillates and even jumps abruptly with varying magnetic flux enclosed 
by the Pb ring, whereas no oscillation was found at positions C and D (data are shown in 
the supplementary materials). 
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FIG. 1: (color online) Illustration of the rf-SQUID used in this experiment. It contains a super¬ 
conducting ring with a gap on the surface of a 3D topological insulator (TI). The proximity-effect- 
induced superconductivity on the TI surface at the gap mediates the Josephson coupling. The 
dots marked by A, B, C, D and E indicate the places where nano-electrodes are attached to the 
TI surface for contact resistance measurement. 

Shown in Fig. 2 are typical results obtained on one of the devices. Figure 2a shows the 
SEM image of the device. The inner and outer diameter of the Pb ring are d, tn = 7.8 /mi 
and d out = 10.2 /m i, respectively. The size of the Josephson junction is 2.2 /tmx200 mu . 
Two Pd electrodes were fabricated and connected to the Bi 2 Te 3 surface at positions A and 
B about 100 nm away from the two ends of the junction, through two windows of 600 nm in 
diameter each on over-exposed PMMA mask. Another Pd electrode was similarly connected 
to the center of the ring (position C). 

Figures 2b and 2c are 2D plots of the contact resistance dV/dl & measured at positions A 
and B, respectively, as functions of magnetic field and bias current. The dV/dlb oscillates at 
a period of A B = 0.32 ± 0.01 G, corresponding to an effective area of S e ff = (fro/ AB = 62.5 
/an 2 (where 0o = h/ 2e is the flux quantum, h is the Planck’s constant, and e the electron 
charge). This estimated area is in good agreement with the geometric area of the ring after 
considering flux compression: -Kd in d out /A = 62.45 /nil 2 . 

The most peculiar feature in the 2D plots is that the low-resistance area has a semilunar 
shape within many oscillation periods, as if an originally 20o-period (i.e., 47T-period) oscil¬ 
lation is truncated when the flux in the ring reaches odd multiples of half flux quantum. In 
addition, the semilunar shapes measured at positions A and B take opposite orientations. 
Within the 0 th envelopes of the Fraunhofer-like pattern (\B\ < 7.7 G), the ones measured 
at position A are centered toward the origin, whereas the ones measured at position B are 
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FIG. 2: (color online) Contact resistance dV/dli between Pd electrodes and Bi 2 Te 3 , measured on 
one typical device near the inner and outer ends of the Josephson junction at T=10 mK. a SEM 
image of the device and the illustration of a three-terminal configuration for contact resistance 
measurement. The black area is covered with over-exposed PMMA which prevents the normal- 
metal Pd electrodes to contact with Bi 2 Te 3 except through small windows at positions A (colored 
in red), B (colored in blue) and C. b and c 2D plots of dV/dlb as functions of magnetic field and 
bias current, measured at positions A and B, respectively. The arrow in c indicates the magnetic 
field at which the Fraunhofer envelope is expected to reach its first minimum. The color scales are 
in the unit of ft. d and e Line cuts of the data in b and c at h,=0, respectively, f The line cuts 
within the dashed window in d and e are plotted together. 

centered away from the origin. The semilunar shapes reverse their orientation after entered 
into the 1 st envelope of the Fraunhofer-like pattern (\B\ > 7.7 G). We note that the char¬ 
acteristic held of 7.7 G corresponds to an area of 2.2 /mix 1.2 /irn, in agreement with the 
effective area of the junction after taking into account of flux compression (explained in the 
supplementary materials). 

Shown in Figs. 2d and 2e are the line cuts in Figs. 2b and 2c along I^ = 0, respectively. 
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With varying magnetic flux, the two sets of data measured at positions A and B vary in 
a complementarity correlated manner, and swap their status abruptly when the flux in the 
ring reaches every odd multiples of half flux quantum. Then the processes start over again. 
The two sets of data in the dashed window in Figs. 2d and 2e are plotted together in Fig. 
2f, showing that there are two 47r-period trends differed by a phase shift of 2tt (i.e., (j) o). 

It has been shown previously that superconductivity can be induced across the Pb- 
Bi 2 Te 3 interface, spreading along the surface in Bi 2 Te 3 up to a distance of microns at low 
temperature^ 13 22 . The measured contact-resistance oscillations at positions A and B there¬ 
fore reflect the variation of the superconducting gap in Bi 2 Te 3 at the ends of the Pb-Bi 2 Te 3 - 
Pb junction. The latter is further controlled by the EPR of the junction. 

The contact resistance of an S-N interface is in general described by the Blonder-Tinkham- 
Klapwijk (BTK) theory^ 23 } When the interfacial barrier is high such that the electron trans¬ 
port across the interface is in the tunneling limit, the contact resistance mainly reflects the 
tunneling density of states of the superconducting side. When the interfacial barrier is low, 
Andreev reflection becomes a dominant process, resulting in a grossly enhanced conductance 
(dirty case) or resonantly enhanced conductance peak(s) (clean case) within the gap energy. 
In our experiment both the high- and low-interfacial barrier cases appeared. In the latter 
case the low-resistance state in the 2D plot could even have a boundary mimicking a critical 
supercurrent, as shown in Figs. 2b and 2c, as if proximity-effect-induced superconductiv¬ 
ity has developed across the Bi 2 Te 3 -Pd interface. This boundary is related to the excess 
current 2 ^ at the S-N interface due to Andreev reflection. 

For the mechanism of jumping, let us firstly point out that it is not related to the 
jumping of the 27r-period screening supercurrent J Si2vr in the ring — a conventional jump¬ 
ing mechanism which is well known to occur when the SQUID’s screening parameter 
/3 e = 2nLI c /(j) 0 > 1 (where L is the inductance of the ring and I c is the critical super¬ 
current of the junction!' 24 ). In fact, for most of the devices investigated in our experiment, 
unless otherwise mentioned, their estimated /3 e was safety below 1 at the base temperature 
(a list of (3 e of our devices can be found in Table SI of the supplementary materials), and was 
further reduced at elevated temperatures due to the decrease of critical supercurrent. For 
devices with j3 e < 1, it is known that their J s 2?r should not jump at any flux. More specifi¬ 
cally, J s 27r should cross zero smoothly at odd multiples of half flux quantum. Therefore, we 
conclude that the dV/dlj, jumping we observed is not caused by the jumping of conventional 
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27r-period screening supercurrent in the ring. This conclusion is further supported by a 
control experiment performed on a device with j3 e = 8.9 (the details are presented in the 
supplementary materials). 

To test whether the dV/dh jumping is unique to the induced topological superconduc¬ 
tivity on TI, we have performed another control experiment on graphite-based devices, in 
which the induced superconductivity in the junction area is supposed to be topologically 
trivial. The results are presented in the supplementary materials as well. It is known that 
the I st 2 n there in the ring must oscillate in its full strength in order to compensate the flux 
change. However, the contact resistance between Pd and graphite responses very smoothly, 
showing no sign of jumping (see Fig. S7-2 in the supplementary materials). It indicates that 
the jumping in dV/dl^ is a unique feature for the devices constructed on Bi 2 Te 3 . 

The abrupt jumping at odd multiples of half flux quantum reflects the existence of a 
27r-period but fully skewed current-phase relation (CPR): I s = I c sin(<^/2) for ip G [—7T,7r]. 
Theoretically, such a CPR could occur in any Josephson junctions made of conventional 
materials, as long as the quasiparticle transport in the junction is fully transparent (i.e., 
with a transmission coefficient D = 1) so that the minigap caused by Andreev refactions 
has a 27r-period form of 24 : A oc |cos(<^/2)|. In practical case, however, the jump will be 
rounded up because D will always be reduced from 1 at finite temperatures and/or due 
to the existence of imperfections. For example, even for a Josephson junction using single 
atomic layer of graphene as the barrier and reached a transmission of as high as D = 0.99, 
its jumping in CPR is already rounded up significantly 25 . For our junctions D is much 
lower. Previous study on interfacial conductances shows that the barrier strength Z of 
Bi 2 Se 3 -Sn interface fabricated with the same technique is around 0.6. It corresponds to 
a transmission coefficient of D = 1/(1 + Z 2 ) = 0.74 - far from being fully transparent. 
Therefore, the appearance of fully skewed and truncated patterns in dV/dl not only at the 
base temperature but also at elevated temperatures (e.g., see Figs. S6-1 and S6-2 in the 
supplementary materials), is rather unusual. 

A fully skewed CPR could arise from a topologically protected mechanism involving two 
branches of EPRs. According to the theories 58 ®^, for Josephson junctions constructed on 
TI surface, there are two branches of 47r-period EPRs when the number of quasiparticle in 
the junction is a good quantum number. Depending on odd or even number of quasipar¬ 
ticles in the system, the two branches belong to two different macroscopic quantum states 
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with odd or even parity. If the system traces the lowest-energy branch at their crossing 
points, then a 27r-period and fully skewed CPR will be yielded. Different from the triv¬ 
ial mechanism containing only one branch of EPR as discussed in the previous paragraph, 
branch-switching between the two EPRs changes the parity of the system, so that it must be 
accompanied with adding or removing one quasiparticle from the system known as quasipar¬ 
ticle poisoning 20 2 ' 28 . Besides the yielded fully skewed CPR, the happening of quasiparticle 
poisoning will also lead to a n phase shift, which can be regarded as a smoking-gun evidence 
whether parity-switching really happens or not in the system. For our device, the number 
of quasiparticles is unfixed because the junction is connected to the quasiparticle bath of 
the surroundings. Therefore, quasiparticle poisoning will happen unavoidably when the two 
branches cross with each other at odd multiples of half flux quantum, if the two-branch 
mechanism really plays a role. 

In the following let us present a detailed interpretation for the boundaries and the enve¬ 
lope of the semilunar shapes shown in Figs. 2b and 2c, as well as for the line cuts shown in 
Figs. 2d, 2e and 2f. It is known that the boundary of the low-resistance state represents a 
characteristic excess current I e caused by Andreev reflection at the BBTe 3 -Pd interface! 23 1 24 1 : 
I e oc (A/eR) tanh(V r /2/csT), where R is the normal-state resistance of the interfacial junc¬ 
tion, and A is the local minigap in Bi 2 Te 3 . For the surface helical electrons with protected 
full transparency, their Andreev reflections in the Pb-Bi 2 Te 3 -Pb junction yield a local mini¬ 
gap of the form: A(x) oc | cos(</?(;r)/2)|. And the local phase difference <p(x ) is further 
controlled by the flux 0 in the ring and in the junction via: 

(p(B,x) = 27 ra^~ + 2 nf3 XH B ± 7rint(-^- + j-) (1) 

00 00 0o 2 

where 0 = BS e ff, x is defined from —W/2 (position A) to W/2 (position B), W = 2.2 /mi 
is the width of the junction, and H* = 1.2 /im is the effective length of the junction after 
considering flux compression (explained in the supplementary materials). The first term in 
Eq. (1) represents the phase difference generated by the flux in the ring. The second term 
is the phase difference generated by the flux in the junction. It determines the slope of the 
the curves in Figs. 3b to 3e. The third term represents an additional n phase shift caused 
by quasiparticle poisoning at every odd multiples of half flux quantum. 

The calculated A at positions A and B are shown in Figs. 3f and 3g, respectively. The 
results pertinently describe the boundaries of the low-resistance state shown in Figs. 2b 
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FIG. 3: (color online) a Formation of a 1/2 fractional mode of Cooper pairs like a Mobius strip on 
the ring, as reflected by the half phase-to-flux ratio ofa = l/2. btoe Distribution and evolution 
with flux of the local phase difference in the junction. The left/right edge of each plot corresponds 
to position A/B of the junction as defined in Fig. 1. The two lines in each plot represent the states 
before and after quasiparticle poisoning. The red dots represent the places where the minigap 
is closed and Majorana fermion states are believed to occur 29 . f and g Calculated minigap at 
positions A and B, respectively, with lineshapes mimicking the boundaries of the low-resistance 
state of dV/dlb shown in Figs. 2b and 2c. h Deduced dV/dlb at /& = 0, mimicking the measured 
line cuts shown in Figs. 2d, 2e and 2f. 

and 2c, with no fitting parameter except for taking a = 1/2 and /3 = 1. A more detailed 
fitting can be found in the supplementary materials. Furthermore, by tentatively taking 
dV/dlb (x e~ A / fcsT (where ks is the Boltzmann constant), the dV/dlb at lb = 0 can also 
be mapped (Fig. 3h). The results mimic the line cuts shown in Figs. 2d, 2e and 2f. In 
particular, the curves within the green window in Fig. 3h correspond to the ones shown in 
Fig. 2f. We note that a rigorous mapping between dV/dlb and A would require numerical 
calculations within the BTK theory®. 
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For ordinary Cooper pairs the phase-to-flux ratio must be a = /3 = 1. The experimental 
finding of a = 1/2 reflects the formation of 1/2 fractional modes of Cooper pairs on the 
ring as represented by the Mobius strip in Fig. 3a. Since the phase generated by the 
flux can in general be expressed a^® (p = 27t0/0q> the half phase-to-flux ratio corresponds 
to an effective flux quantum of </>q = h/e* = 20o ; where e* = e is the effective charge of 
fractionalized Cooper pairs that passes through the junction, mimicking in a superconducting 
Kitaev chain 7 . It is the appearance of a = 1/2, not the full-transparency-caused 1/2 factor 
in the minigap formula, that provides the evidence for the existence of the long-sought 
47r periodicity, despite that the measured periods are 27r-periodic after being truncated by 
quasiparticle poisoning. 

Contrary to the global mode on the ring which passes through the TI surface once to form 
the 1/2 fractional mode, the local mode surrounding the Josephson vortex in the junction 
passes through the TI surface twice. We believe that this guarantees the formation of integer 
modes in the junction, so that the phase-to-flux ratio for the second term in Eq. (1) is /3 — 1. 
More theoretical study would be needed to discuss this experimental finding. 

According to Potter and FiP^, the red dots in Figs. 3c to 3e represent the positions where 
the minigap is closed so that Majoranas are expected. The positions of Majoranas can be 
manipulated by varying the flux, and swapped by quasiparticle poisoning when the two dvr- 
period EPRs cross with each other at odd multiples of half flux quantum. Figures 3c and 
3e further explain that, when the flux in the junction area reaches half flux quantum (i.e., 
~ 11.50o bi Fig. 3), the Majoranas move to the ends of the junction, so that the contact 
resistance at position A/B experiences full gap-closing, together with the sharpest jumping 
before and after position swapping. 

To conclude, we have observed clear evidence for the existence of 1/2 fractional modes and 
the occurrence of quasiparticle poisoning. These are necessary signatures accompanied with 
the formation of Majorana fermion states. Our study demonstrates that superconducting 
devices constructed on the surface of 3D TIs provide a promising platform for braiding 
Majorana fermions in two dimensions and building topological quantum computers in the 
future. 


* These authors contribute equally to this work. 
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1. Contact resistance between Pd and Bi 2 Te 3 at other positions of the rf-SQUID 

We have fabricated several devices which contain Pd electrodes not only at positions A and 
B, but also at positions C and D defined in Fig. 1 of the main manuscript. Figure Sla shows the 
SEM image of one of such devices (device #S1). The 2D plot of contact resistance dV/dl b 
measured at positions A, C, and D are shown in Figs. SI b, c, and d, respectively. Only at positions 
near the line Josephson junction, e.g., at position A, dV/dl b oscillations were observed with 
varying magnetic flux in the Pb ring. No resistance oscillations were found at positions away from 
the line Josephson junction, namely positions C and D, no matter their positions are close to the 
superconducting ring or not. The results indicate that the oscillations we observed are only 
related to the status of the Josephson junction in the rf-SQUID loop. 

Although there is no oscillation, the dV/dl b data measured at position D still show a 
low-resistance state within a characteristic gap energy, as shown in Fig. Sid. This low-resistance 
state is caused by the induced superconductivity in Bi 2 Te 3 in the vicinity of the Pb ring. In the 
main manuscript, we attribute the low-resistance state to Andreev reflection and/or proximity 
superconductivity at the secondary interface (i.e., between Pd and Bi 2 Te 3 ). 

In addition to the low-resistance state, there is also a zero-bias resistance peak (ZBRP) in Fig. 
Sid. This ZBRP is a common feature of the secondary interface, regardless there are oscillations 
or not — the ZBRP superimposes with the dV/dl b oscillations at positions A and B when the 
contact resistance is relatively high (see Fig. S2-lb, also Figs. S2-2a and S2-2b). This ZBRP at the 
secondary interface is directly related to the zero-bias conductance peak (ZBCP). The latter is 
reproducibly observed at the primary interface [1], It is presumably related to the gap-closing in 
Bi 2 Te 3 at the Fermi energy. 
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FIG. SI | (a) SEM image of device #S1. (b), (c) and (d) The dV/dl b measured at positions A, C 
and D, respectively, as functions of magnetic field and bias current at T=10 mK. 
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2. Data obtained on Pd electrodes with relatively high contact resistance 


As mentioned in the main manuscript, for the data taken by electrodes with relatively low 
contact resistance, there is a clear boundary in the 2D plots of dV/dl b mimicking a critical 
supercurrent for the low-resistance state, indicating that proximity superconductivity has 
developed across the S-N interface. 

We have also investigated a number of devices with relatively higher contact resistance 
between Bi 2 Te 3 and Pd, but still in the regime of conductance enhancement within the gap. We 
found that their low-resistance state no longer has a clear semilunar boundary anymore. The 
absence of such a boundary indicates that proximity superconductivity is not yet developed 
across the Bi 2 Te 3 -Pd interface, only conductance enhancement via Andreev reflection plays a role 
there. Nevertheless, semilunar-like shapes still emerge in the 2D plots. 

In Fig. S2-1 we show the data obtained on device #S2 which has electrodes with both 
relatively high and low contact resistance. The data in Fig. S2-lc, where the low-resistance area 
has a clear boundary, indicate that the S-N interface at position B has become superconducting. 
This argument is further supported by the fact that the characteristic critical supercurrent 
increases with decreasing temperature, as will be shown later in Figs. S6 f, g, h, I, and j. On the 
other hand, for the data shown in Fig. S2-lb, the low-resistance area does not have a clear 
boundary. It indicates that proximity superconductivity is not yet well developed across the S-N 
interface at position A, only conductance enhancement via Andreev reflection occurs. This 
argument is again supported by the temperature variation shown in Figs. S6 a, b, c, e, and f. With 
decreasing temperature, the enhanced conductance peak gets narrower and narrower against 
both bias current and bias voltage, exhibiting a resonant nature. 



FIG. S2-1 | The contact resistance dV/dl b of device #S2 measured at 10 mK. (a) An SEM image of the 
device, (b) 2D plot of dV/dl b measured at position A. The oscillating low-resistance region does not 
have a clear boundary, (c) 2D plot of dV/dl b measured at position B. The low resistance regions have 
a clear boundary. 


Regardless of relatively high or low contact resistance, the jumping of dV/dl b with magnetic 
flux can always be seen at positions A and B, both on the line cuts of the 2D plot and through the 
hyperbolic-like and concentric-like semilunar shapes in the 2D plot. For electrodes with high 
contact resistance, the jumps on the line cuts appear to be even sharper. When the magnetic flux 
is slightly away from zero, the contact resistance periodically oscillates to its normal-state value, 
reflecting that gap-closing happens in Bi 2 Te 3 at odd multiples of half flux quantum. 

In Fig. S2-2 we show more data obtained on device #S3 whose both electrodes at positions 
A and B have relatively high contact resistance. There are sharp jumpings to the normal-state 
resistance (namely gap-closing) on the line cuts. 
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FIG. S2-2 | Oscillatory patterns of dV/dl b obtained on another device (device #S3) with relatively 
high contact resistance. The measurement was performed atT=10 mK. (a) and (b) dV/dl b measured 
at positions A and B, respectively. Although without a clear boundary, hyperbolic-like and 
concentric-like semilunar shapes still emerge, (c) and (d) Line cuts in (a) and (b) at l b = 0, 
respectively. Gap-closing in Bi 2 Te 3 happens at odd multiples of half flux quantum, as that the dV/dl b 
hits the normal-state values (the dashed lines) of the contact resistance at these flux. 


3. Discussions on p e and a list of p e for the devices investigated 

In applied magnetic field, screening supercurrents will be induced in our device, circulating 
along the Pb ring and around the Josephson junction area. When the critical supercurrent of the 
junction is large and/or the inductance of the loop is large, the induced supercurrent in the ring 
can significantly influence the total magnetic flux in the ring, trying to make it quantized. 

The total magnetic flux in the loop is [2]: 

<5 = <5e - (Pe<l>o/2Jt) sin(2jt<J>/<t> 0 ) (1) 

where <l> e is external flux exerted by the applied magnetic field, (|)o=h/2e is flux quanta, 
p e =2xrU c /(po is the SQUID screening parameter, L is the inductance of the ring, and l c is the critical 
supercurrent of the junction in the superconducting loop. 



FIG. S3 | The total magnetic flux €> vs. the external magnetic flux <t> e for SQUID loops with 
different SQUID screening parameter p e . (a) p e <l. (b) p e =l. (c) p e >l. For the case (c), flux 
jumping, 2jt-period circulating supercurrent jumping, and hysteresis between opposite 
directions of field sweeping occur. 


4 





























The dependence of the total flux <J> on the external flux ® e shows two different kinds of 
behavior according to the value of p e - For p e <l, < 1 ) is single valued and increases monotonically 
with increasing ® e (Fig. S3 a and b). For p e >l, ® is multivalued on ® e , only the segments of the 
curve with positive slope are traced in field sweeping (Fig. S3c). There is hysteresis between 
opposite directions of field sweeping. The total flux of the ring and the circulating supercurrent in 
the ring jump at positions determined by the specific value of p e , not necessary at odd multiples 
of half flux quanta. 

To determine the p e of our rf-SQUID, we need to know the critical supercurrent of the 
Josephson junction in the SQUID loop which cannot be directly measured because the junction is 
short-circuited with the superconducting loop. We thus have to fabricate a single Josephson 
junction of the same size and on the same flake of Bi 2 Te 3 or graphite, then using the critical 
supercurrent of the single Josephson junction to estimate the p e of the SQUID. Some of the 
results are summarized in Table SI. 


Table SI | The (3 e of some of the devices we investigated 


Device number 

Device 
presented in 
the main 
manuscript 

#S2 

#S3 

#S4 

#S5 

Geometric 
parameters of 
the Pb rings and 
the junctions 

D in =7.8 pm 
D out =10.2 pm 
W=2.2 pm 

D in =5.7 pm 
D out =8 pm 
W=2.2 pm 

D in =5.7 pm 
D out =8.1 pm 
W=3.2 pm 

D in =4 pm 
D out =5.8 pm 
W=2.2 pm 

D in =18 pm 
D ou ,=20 pm 
W=5.2 pm 

Inductance L 

9.9 pH 

5.7 pH 

5.8 pH 

2.9 pH 

44.8 pH 

Estimated Ic of 
the junction 

~ 27 pA 

~ 27 pA 

~ 39 pA 

~ 27 pA 

~ 63 pA 

Estimated (3 e 

0.84(5)10 mK 

0.48(6)10 mK 

0.71(6)10 mK 

0.25(6)20 mK 

8.9(6)20 mK 


Notes: D in and D out are the inner and outer diameters of the Pb ring, respectively. 

W is the width of the Josephson junction. The length FI of the junctions is the same for all 
devices, being around 200-300 nm. 


For most of the devices we investigated, p e is smaller than one, so that the devices are in the 
non-hysteretic regime. For example, the data taken on device #S4 are shown in Fig. S4c. No 
hysteresis is seen between the black and red lines taken in opposite directions of field sweeping. 
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4. More data obtained on devices with p e <l 


Shown in Fig. S4 are the data obtained on another device #S4 whose estimated SQUID 
screening parameter is p e =0.25. It has a slightly smaller diameter and also a slightly higher 
contact resistance than that of the device shown in Fig. 2 of the main manuscript. The 
temperature and field dependencies of the data clearly indicate that it is the resistance dip (i.e., 
the conductance peak) near zero bias voltage that oscillates with the magnetic flux in the ring. 
The energy scale (width) of this dip that oscillating is around 0.2 mV, an order of magnitude larger 
than that of the ZBRP structure shown in Fig. Sid, indicating that this dip is an 
Andreev-reflection-induced gross conductance enhancement within the gap of superconducting 
Bi 2 Te 3 . Proximity superconductivity is not yet developed across the S-N interface, because the 
low-resistance state in Fig. S4b does not have a clear boundary. The amplitude of resistance 
oscillation in Fig. S4c corresponds to an conductance enhancement of 0.45 e 2 /h. Figures S4 e and 
f show that, with decreasing temperature, the width of the dip gets narrower, and the depth of 
the dip gets deeper, showing a resonant nature. This indicates again that the low-resistance state 
is caused by Andreev reflection, not the appearance of a supercurrent. 


S o 



B(G) 






B(G) 


T(mK) 


FIG. S4 | The results obtained on another device #S4 with p e =0.25. (a) An scanning electron 
microscopy image of the device. The single Josephson junction identical to that on the ring was 
used for estimating the critical supercurrent of the junction, (b) 2D plot of dV/dl b measured at 
position B (near the end of the junction outside the ring) at T=30 mK. (c) dV/dl b at line cut l b =0, 
measured in two opposite field-sweeping directions (lines in black and red, respectively), 
showing no hysteresis. Also, the curves periodically hit the normal-state value represented by 
the dashed line, indicating gap-closing in Bi 2 Te 3 at these flux, (d) The bias voltage dependence of 
dV/dl b measured at position B at T=30 mK and in different magnetic fields, (e) The temperature 
evolution of the bias voltage dependence of dV/dl b measured at position B at zero field, (f) The 
temperature dependencies of the depth (black squares) and the full width at half depth (red 
circles) of the dip. The ac excitation current was 20 nA except in (c) where it was 2 nA. 
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5. Control experiments (I): data obtained on a devices with p e >l 


To demonstrate that our estimation of p e is reliable, we have performed a control 
experiment on a device with p e > 1. In this case, l s , 27 t is expected to jump via the conventional 
mechanism, at the places which is not necessary located at odd multiples of half flux quantum 
but rather depending on the specific value of p e . And hysteresis in jumpings is expected if the 
field sweeping direction is reversed. 

In Fig. S5 we show the data obtained on device #S5 which has a large loop area and a long 
junction width (5.2 pm). Its screening parameter is estimated to be p e =8.9. Obvious hysteresis 
happens when the field sweeping direction is reversed. 

We must point out that in the conventional mechanism there should have no oscillation and 
jumping in the hysteresis region. The jumping observed here in this region is presumably still 
caused by quasiparticle poisoning, the accumulation of which could even prevent the happening 
of the conventional jumping. Nevertheless, the position of jumping now is influenced by the 
screening supercurrent via the conventional mechanism, so that the places of jumping for this 
device are not located at every odd multiples of half flux quanta, but near the integer multiples of 
flux quanta. This is understandable. If the conventional screening supercurrent l S/2jI is large 
enough in large-p e device, and is able to modify the total magnetic flux in the ring, then the 
places of jumping caused by the anomalous l Sj4ll supercurrent will be shifted accordingly. 

On one hand, the results tell us that the places of jumping in dV/dl b can be influenced and 
shifted away from odd multiples of half flux quantum in the p e > 1 case. On the other hand, the 
results confirm that the rest devices we investigated are all in the p e < 1 regime, for that their 
jumpings always take place precisely at odd multiples of half flux quantum, showing no hysteresis 
in bi-directional field sweeping (e.g., see Fig. S4c). Therefore, we can conclude that the dV/dl b 
jumping we observed is not caused by the jumping of the conventional 2n-period screening 
supercurrent in the ring. 



FIG. S5 | dV/dl b curve at position A of device #S5 (P e =8.9) as a function of magnetic field, 
measured in opposite field sweeping directions (illustrated by the arrows). T=20 mK. 


7 






















6. The temperature dependence of the contact resistance oscillation 


To further support our argument that the observed jumping in contact resistance is not 
caused by the jumping of the conventional 2xc-period supercurrent l s , 2lI in the ring in the p e >l 
case, we have performed measurements on several devices not only at the base temperature, but 
also at elevated temperatures, to further reduce the l c and thus the p e . 

Figure S6-1 shows the contact resistance at positions A and B measured on device #S2 
(Pe=0.54 at 10 mK) at several different temperatures. Although p e gets smaller and smaller with 
increasing temperature, the jumping phenomenon persists and remains to be sharp (Fig. S6-2). 
Moreover, the places of jumping keep locked at odd multiples of half flux quanta for the data 
taken at different temperatures. The results convincingly rule out the possibility that the jumping 
in contact resistance is caused by the jumping of the conventional 2n-period supercurrent l s _ 2lt in 
the p e >l case. 

On the other hand, if the occurrence of jumping corresponds to a fully skewed 
current-phase relation of our devices (i.e., a fully transparent barrier), then it is quite surprised 
that the jumping keeps sharp at elevated temperatures to T/T c =0.4 - 0.6, as shown in Fig. S6-2. 
The jumping caused by this mechanism usually rounds up with increasing temperature and/or 
disorders. The results provide strong evidence that the jumping has a non-trivial mechanism. 


T=10mK 

Kill 


-0.2 


b 50 mK 

: til 


-0.2- 

0.2- 

< 

si 0.0 

-0.2 


d 150 mK 

11 Itll 


-0.2 


e 200 mK 

1(1111 


-0.2 


100 mK 

III 


dV/d l b 
(O) 

H 295.0 
I 286.3 

m 2776 

268.9 
— 260.2 



B (G) 


dV/d l b 

(Q) 

1 68.0 

65.7 
63.4 
61.2 
— 58.9 


FIG. S6-1 | Temperature 

dependence of the contact 
resistance dV/dl b measured 
at position A (a, b, c, d, e), 
and at position B (f, g, h, i, j) 
on device #S2. 












FIG. S6-2 | (a) Line cuts at l b =0 in Figs. S6-1 (a), (b), (c), (d) and (e), respectively, (b) Line 
cuts at l b =0 in Figs. S6-1 (f), (g), (h), (i) and (j), respectively. The jumpings remain to be 
sharp at elevated temperatures of T/T c = 0.4 - 0.6, where T c = 0.3 - 0.5 K is the critical 
temperature of proximity-induced superconductivity in Bi 2 Te 3 . 
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7. Control experiments (II): data obtained on graphite-based devices 


We have also performed control experiments on graphite-based devices. The results are 
shown in Figs. S7-1 and S7-2. 

Unlike the electron system on the surface of a 3D Tl which contains only one type of helical 
electrons, the electrons in graphene are four-fold degenerated. Therefore, Graphite is known as a 
topologically-trivial material, with negligible spin-orbit coupling. It is an ideal candidate to be 
used for performing comparative measurements. 

The superconducting proximity effect between Pb and graphite appears to be relatively 
weak. We therefore have to use Sn to replace Pb. Due to various technique issues, it took us 
several months to succeed. 

In Fig. S7-1 we show the results obtained from one of the graphite-based devices. The 
contact resistance measurement by Pd electrode at position B revealed a fully developed gap at 
low temperatures (i.e., with saturated width and amplitude, see Fig. S7-lb), reflecting that the 
graphite beneath the electrode has become superconducting, owing to the proximity effect from 
the Sn rf-SQUID. In such a device, the l s2n in the ring must oscillate in its full strength trying to 
compensate the change of magnetic flux in the ring. Flowever, no noticeable influence on the 
contact resistance (hence, the gap) was observed during this process, only a Fraunhofer-like 
pattern of the Josephson junction was seen (Fig. S7-lc). 

The Fraunhofer-like pattern shown in Fig. S7-lc is for a high resistance state. It reflects that 
the gap beneath the Pd electrode in graphite is modulated by the flux in the Josephson junction. 
A similar Fraunhofer-like pattern for a zero-resistance state was observed on the single Josephson 
junction in the right side of Fig. S7-la (data not shown). 



-40 -20 0 20 40 -20 0 20 

V (nV) B (G) 

FIG. S7-1 | A comparative experiment carried out on Sn-graphite device, (a) An scanning 
electron microscopy image of the device, (b) The bias voltage dependence of the contact 
resistance dV/dl b measured at position B at several different temperatures and in zero 
magnetic field, (c) 2D plot of dV/dl b measured at position B at T=10 mK, demonstrating a 
Fraunhofer-pattern-like variation of gap with magnetic flux in the Josephson junction area. 


One possibility of observing no contact resistance oscillation but only a Fraunhofer-like 
pattern is that the probing current l b (which was applied between the Pd electrode and the Sn 
ring) was large compared with the critical supercurrent of the junction, so that it disturbed the 
current flowing in the junction. This could happen because the parameters such as l c for devices 
on graphite were not well controlled so far — the junction on the ring might have a small l c 
compared with the applied l b . We therefore investigated more devices. 
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FIG. S7-2 | Results obtained on another Sn-graphite device, (a) An scanning electron 
microscopy image of the device, (b) 2D plot of dV/dl b measured at position B at T=10 mK. Note 
that the range of the color scale is only 1% of the total resistance, much smaller than the ones 
used for other devices. 

In Fig. S7-2 we show the results obtained on another graphite-based device. An oscillatory 
pattern can be marginally resolved in a period corresponding to the ring area. However, the 
amplitude of dV/dl b oscillation is less than 0.3%, about two orders of magnitude weaker than 
those observed on Tl-based devices. And no sign of jumping can be resolved. 

Overall, our results indicate that the superconductivity based on topologically-trivial bulk 
states gives no contribution to the dV/dl b jumping. Thus, the jumping observed on Tl-based 
rf-SQUIDs is presumably a phenomenon particularly related to the superconducting surface of 
Bi 2 Te 3 . 

In the following we give an explanation for the observed result shown in Fig. S7-2b. Similar 
to the analysis in the main manuscript, we assume that what the Pd electrode probed is the 
minigap in graphite at position B near the Josephson junction of the rf-SQUID. And we assume 
that the quasiparticle transport within the junction area is nearly ballistic so that the minigap can 
still be expressed as Aoc | cos(cp/2) |. Because in graphite the electron states are trivial, the 
phase-to-flux ratio of the ring is a=l. Therefore, cp can be expressed as: 
cp(B, x) = 27Ccj)/<|)o + 2xcxHB/(j)o 

where B is the magnetic field, x is defined from -W/2 (position A) to W/2 (position B), W is the 
width and H is the length of the junction. 

Figure S7-3 shows the calculated pattern of minigap oscillation at position B (at x=W/2), 
which is in good agreement with the experimental result shown in Fig. S7-2c as long as the 
characteristic current in Fig. S7-2c is proportional to the minigap. 



FIG. S7-3 | Calculated pattern of minigap oscillation at position B (x=W/2) 
in graphite-based rf-SQUIDs. 

Because electron transport in real graphite devices is not fully transparent, the minigap will 
never be fully closed, which accounts for the weak oscillation amplitude of dV/dl b . 
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8. Josephson energy profiles in the presence of both 2n-period and 4n-period modes, 
quasiparticle poisoning and supercurrent reversing 

In many phase sensitive experiments, the trivial 2xc-period CPR, which would mainly arise 
from the bulk states, gives considerable contribution. In the presence of both 2n-period and 
47t-period modes in a rf-SQUID, the total Josephson energy is the sum of these modes. 

For the conventional 27i-period mode, if assuming that the magnetic energy of the ring 
inductance is small (which is true in our case because of small critical supercurrent and ring area), 
then the total energy is: 

Ej,2ti °c (VW) f [l-cos( 27 i(()/(t)o+ 27 iBHx/(t)o)]dx (for x=-W/2 to W/2) 

= -A sin{7i(|)'/(t)o)/(7t(|)'/(|)o)cos(27i(|)/(t)o) 

where ((> is the magnetic flux in the ring, and (|)'=BHW is the magnetic flux in the junction, B is the 
magnetic field, H is the length and W is the width of the junction. 

Similarly, for the two 47t-period modes (a=J/2) with a 2n phase shift, their energy can be 
calculated by using the phase expressed in Eq.(l) of the main manuscript. The results are: 

Ej, 4n = + e sin(7n|)'/(|)o)/(7i(j)'/(|)o)cos(7i(j)/(|)o) 
where the plus/minus sign represents the odd/even parity branch. 

Thus the total energy of the device can be written as: 

Ej = sin(7i<|) , /(t) 0 )/(7i<t) , /(t> 0 ) [—A cos(27T(j)/(J) 0 ) i B cos(7i<))M>o)] 



- 1 - 1 - i - 1 - i - 1 - i - 

- 2-10 1 2 


FIG. S8 | Energy profiles of a rf-SQUID in the presence of both 2n- (black) and 47i-period 
modes (dashed red and blue) in the vicinity of The total energy is represented by the 
solid red and blue curves. The arrows indicate where branch-crossing happens. 

In Fig. S8, E lp27t is plotted in black line, and E iM is plotted in dashed red and dashed blue lines 
for even and odd parity modes, respectively. The total energy of the system is plotted in solid red 
and solid blue lines, assuming that 6=0.54. 

It can be seen that the place where the two branches cross with each other is always located 
at odd multiples of half flux quantum for whatever ratio of B/A. There the slope of E Ji2ji is zero, so 
that the physics is dominated by E Jj4lt . 

The crossing of the two branches allows the happening of quasiparticle poisoning, which will 
be right at odd multiples of half flux quantum. 

Given the fact that the supercurrent is proportional to the first-order derivative of the 
energy, quasiparticle-poisoning-induced branch-switching causes no effect on the 2n-period 
component of supercurrent, but reverses the 4n-period supercurrent. 
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9. Evolution of the phase and the minigap with magnetic flux, more details 


In the following we give more plots to illustrate the evolution of the phase (thus the gap) 
with magnetic flux in the ring and in the junction area, in the presence of the fractional modes. 

We note that a flux of 4 > 0 fl corresponds to a phase of n if a=l, or to ji/ 2 if a=l/2. In the 
absence of the 4n-period modes, the 2n-period mode should undergo gap-closing at odd 
multiples of (j)o fl. In the presence of the 47i-period modes, however, the gap actually does not 
close at odd multiples of c() 0 /2 when the magnetic flux in the junction is less than 4> 0 /2 [i.e., in the 
regime of | B| <3.8 Gauss in Figs. 2(B) and 2(C)], Nevertheless, phase jumping still occurs in this 
regime (see Fig. 2 of the main manuscript), indicating the occurrence of quasiparticle poisoning 
due to the existence of two 4n-period branches of energy-flux relations and their degeneracy 
right at odd multiples of (j) 0 /2, as discussed in the previous section. 

As to the positions of the Majoranas, we anticipate that Majoranas might be delocalized 
inside the junction in the regime of | B | <3.8 Gauss [e.g., for the case illustrated in Fig. S9-l(C)j if 
they exist. In this regime the junction could be regarded as a short junction, as treated by Wieder, 
Zhang and Kane [3], Although the minigap in this regime is not closed at odd multiples of half flux 
quantum because of the half phase-to-flux ratio, the two 47r-period EPRs of the system still cross 
with each other (as discussed in Section 8 and shown in Fig. S8), yielding superpositions of even 
and odd-parity branches just like the superposition of occupied and unoccupied quasiparticle 
states. Quasiparticle poisoning could still happen at the degenerate points, as has been observed. 

In the regime of | B | >3.8 Gauss, the increased flux in junction tilts the phase-position curves, 
letting them to reach odd multiples of n where the minigap is closed [see Figs. S9-2(C)d, S9-3(C), 
S9-4(C)j, so that Majoranas as superpositions of occupied and unoccupied states become 
restricted in spaces in the junction, as discussed by Potter and Fu [4]. The positions of the 
Majoranas are marked by the red dots in Figs. 3 (B) to (E) in the main manuscript, and also in the 



FIG. S9-1 | (A) and (B) Gap energy at positions A and B. (C) a to e Phase distribution in the 
junction, at ring flux indicated by the gray lines, cp may differ from the real value by integer 
multiples of 2ji depending on the details of quasiparticle poisoning. 
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FIG. S9-2 | (A) and (B) Gap energy at positions A and B. (C) a to e Phase distribution in the 
junction, at ring flux indicated by the gray lines, cp may differ from the real value by integer 
multiples of 2n depending on the details of quasiparticle poisoning. In d, Majoranas start to be 
restricted in the junction whose positions are marked by the red dots. 
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FIG. S9-3 | (A) and (B) Gap energy at positions A and B. (C) a to e Phase distribution in the 
junction, at ring flux indicated by the gray lines, cp may differ from the real value by integer 
multiples of 2 ji depending on the details of quasiparticle poisoning. 


14 


















































































































UnA) 



FIG. S9-4 | (A) and (B) Gap energy at positions A and B. (C) a to e Phase distribution in the 
junction, at ring flux indicated by the gray lines, cp may differ from the real value by integer 
multiples of 2n depending on the details of quasiparticle poisoning. 


10. Detailed comparison between the data and the model 

The minigap shown in Fig. 3 of the main manuscript is obtained with no fitting parameters 
except for letting a=l/2 and (3=1. For a more detailed comparison between the data and the 
model, we need to put in some fitting parameters. 

In the following, let us fit the boundary of the low-resistance state in the 2D plot of dV/dl b 
with our model. We assume that the boundary corresponds to a total characteristic current l c ' 
which contains the contributions of an excess current l e (l e ccA), and an empirical leakage current 
(being a constant, presumably leaking into the bulk superconducting state): 

IcV Ico' = aA/Ao + b 

where l c0 ' is the maximum value of total characteristic current, A 0 is the maximum value of 
minigap, a and b are fitting parameters. 

By taking a/b = 3 and using the minigap presented in Fig. 3g of the main manuscript, we 
obtained the yellow curve in Fig. S10. The agreement between the data and the model seems 
remarkably well. 

In the main manuscript, nevertheless, we prefer just to present the A vs. flux curves in Fig. 3 
for simplicity. 

180 


150 
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— 90 

FIG. S10 | Comparison between the model (the yellow line) and the boundary of the low-resistance 
state shown in Fig. 2c of the main manuscript. 
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11. On the position uncertainty of the Majoranas and the flux uncertainty at which the jumps 
happen 


(1) On the position uncertainty of the Majoranas 

In Wieder et al.'s theory [3] the flux in the junction area is not considered (i.e., in the 
low-field limit, or in the small junction approximation). In this case the Majorana is presumably 
delocalized inside of the junction area. 

In Potter and Fu's theory [4], on the contrary, the flux in the junction area is considered (i.e., 
in the large-field limit, or in the large junction approximation). The Majoranas are believed to 
exist at places where the local phase difference crosses n such that the minigap is closed. 
However, due to thermal smearing, these places have an uncertainty. From the data of the device 
shown in Fig. 2 of the main manuscript we know that its minigap oscillate between 0 and ~28 pV. 
Assume a thermal energy of AE T =2 pV at an electron temperature of 20 mK, the phase 
uncertainty 25 can be estimated from the relation AE T = 28pV| cos[(jt-5)/2] |. We got 5 =0.0467t. 
Such a phase uncertainty corresponds to a length scale uncertainty of the Majorana Ax = 25W = 2 
X 0.046 X 2.2 pm = 0.2 pm when the phase difference between A and B is around ji (e.g., for the 
case of Fig. 3c of the manuscript). At higher flux the position uncertainty becomes smaller, 
because of the larger slope of the phase-position curves. 

(2) On the flux uncertainty at which the jumps take place 

As we have pointed out in Section 8 of the supplementary materials, quasiparticle poisoning 
will happen at odd multiples of half flux quantum where the even- and odd-parity branches cross 
with each other. Since the energy scale of the two branches is not but the minigap, the phase 
uncertainty 25 =2 X0.04671=0.092n also defines a flux uncertainty window of ~ 0.1())o within 
which quasiparticle poisoning (thus jumping) can take place. From our data, however, the 
uncertainty window for the happening of quasiparticle poisoning appears to be smaller (being a 
few tenth of 0.1(|) 0 ), which is understandable because in our measurement the flux is ramped 
slowly along one direction. 


12. Why is the 47t-period energy-phase relation the dominant signal measured? The role of the 
minigap of the surface states versus the superconducting gap of the bulk states 

In our Bi 2 Te 3 samples there are both helical surface states and bulk states. 

The helical electrons on Bi 2 Te 3 surface have relatively long mean-free-path, and presumably 
with protected transmission at the Pb-Bi 2 Te 3 interfaces, so that they undergo multiple Andreev 
reflections between the two Pb electrodes of the Josephson junction. The yielded local Andreev 
bound states further form quasi ID structures along the width direction of the junction. The 
lowest Andreev bound state defines the size of the minigap. And the minigap is a function of 
position and magnetic flux in the junction. 

The bulk states also become superconducting due to the proximity effect at low 
temperatures, but forming no Andreev bound states at finite energies in the energy window of 
our experiment, as reflected by the experimental data, presumably because of their trivial nature. 
They would remain to be superconducting while the minigap of the surface state is tuned to be 
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closed by the magnetic flux in the junction. 

The contact resistance measurement in our experiment is most sensitive to the surface state 
of the Tl. It mainly probes the oscillation and closing of the minigap of the lowest Andreev bound 
state on the Tl surface, despite that the bulk superconducting gap and the related 2n 
supercurrent might response differently. This probably explains why a fully skewed EPR is 
observed in our contact resistance measurement, but a fully skewed CPR was not observed in 
previous interference/Fraunhofer pattern measurements for critical supercurrent, of which the 
bulk superconducting state also contributes. 


13. On the three-terminal measurement configuration for contact resistance measurement 


In Fig. S13 we re-plot the rf-SQUID which shows the positions of electrodes A, B, C, and E. 
We have also sketched on top of the figure the current flowing paths (the red lines) and voltage 
measurement path (the green lines) in a three-terminal measurement configuration. 



FIG. S13| Left panel: Current flowing and voltage drop in our three-terminal measurement 
configuration. The red lines represent the flowing of normal current, and the blue lines 
represent the flowing of supercurrent. The current injected from electrode A tends to flow 
into the superconducting ring, redistributing on the ring, then flowing from the ring to 
electrode C. The green lines represent the voltage measurement loop. Right panel: the 
equivalent circuit in a cross-section view. 


What is the signal we measured in the three-terminal configuration? 

Our general picture is: the flux determines the local phase difference in the Pb-Bi 2 Te 3 -Pb 
junction, and the local phase difference determines the local supercurrent density (via the 
Josephson equation) as well as the local minigap in the junction. The local minigap further 
determines the measured contact resistance dV/dl b at the Bi 2 Te 3 -Pd interface. 

From the right panel of Fig. S13 it can been seen that along the voltage measurement loop 
there are voltage drops: (1) on R P d fii m ,A which is <30 Q and being a constant; (2) on R P dcontact,A 
which ranges from ~ 100 Q to several kQ in our experiments, and partially oscillating with the 
flux; and (3) on Bi 2 Te 3 surface which is small, with an equivalent non-local resistance of only ~ 1 Q 
(see Y. Pang, et al., arXiv:1503.00838vl). In such case, the dominant oscillatory signal we 
measured in the three-terminal configuration, with oscillation amplitude of ~ 10 Q along the line 
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cut at zero-bias current and up to several tens Q along line cuts at high bias currents, can only 
come from the oscillation of the contact resistance. 

Indeed, in the beginning of the experiment we have checked the contact resistance 
measurement with different combinations of electrodes, such as to measure the contact 
resistance of electrode A with whatever electrodes C, D, E or even the ring as the second and 
third electrodes, the semilunar shape and jumping looked to be identical. 


14. Notes on the measurement currents 

In our experiment, the ac excitation current used to measure the dV/dl b is around 1 nA, 
being four order of magnitude smaller than the l c of the Pb-Bi 2 Te 3 -Pb junction (several 10 pA). 

The ramping range of dc bias current l b was ±0.1 to ±1 pA, about two order of magnitude 
smaller than the l c of the Pb-Bi 2 Te 3 -Pb junction which was a few tens pA for the device shown in 
Fig. 2 of the manuscript. Therefore, the ac and dc currents used in our measurement would not 
influence the status of the Pb-Bi 2 Te 3 -Pb junctions. 

Further increasing l b to the l c of Pb-Bi 2 Te 3 -Pb junction will ultimately influence the status of 
the junction, driving it to the normal state. We believe that this is the reason why we observed a 
Fraunhofer-like pattern of contact resistance on one of the graphite-based devices, where we 
used the superconducting ring as one of the current leads. 

With the ramping of l b , the applied bias current eventually exceeds the local characteristic 
supercurrent (several 0.1 pA) of Bi 2 Te 3 which is determined by the minigap of the junction. In this 
way the amplitude of the local minigap is measured. 


15. Estimation of the effective junction area in the presence of flux compression, stray 
supercurrent and proximity-effect-induced superconductivity on Bi 2 Te 3 surface 

In Fig. S15 we illustrate how the effective junction area of the device shown in Fig. 2 of the 
main manuscript is estimated in the presence of flux compression, stray supercurrent distribution 
and proximity-effect-induced superconductivity on Bi 2 Te 3 surface. Our previous study shows that 
the proximity-effect-induced superconductivity spreads from the Pb electrodes to a distance of 
micron on the surface of Bi 2 Te 3 . (Qu, F. et al., Sci. Rep. 2, 339 (2012); Yang, F. et al., Phys. Rev. B 86, 
134504 (2012)). For this reason, at the two ends of the junction there will be areas distributed 
with stray supercurrents on the surface of Bi 2 Te 3 . Since the strength of superconductivity in these 
stray areas decays with the distance, their edges are not well defined. Adding to the complexity, 
the effective area might also vary with the local supercurrent density. Nevertheless, with our 
experience on the Fraunhofer pattern of single Josephson junctions of the same kind and with 
similar size, we came up with the empirical method and understanding as shown in the right 
panel of Fig. S15 for estimating the approximated effective area of the junctions. The error is less 
than 20%. 
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FIG. S15| Upper left panel: Illustration of flux compression. Lower left panel: The flux within each 
area defined by the white dashed lines in Pb will be compressed along the arrow direction. Right 
panel: Estimation on the effective junction area in the presence of flux compression and stray 
supercurrent distribution on Bi 2 Te 3 surface. The effective junction area for the device shown in 
Fig. 2 of the main manuscript can be approximately represented by the red box, which is 1.2x2.2 
pm 2 . 
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